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Abstract 

Nitrogen vacancy (NV _ ) color centers in diamond are a prime candidate for use in quantum 
information devices, owing to their spin-1 ground state, straightforward optical initialization and 
readout, and long intrinsic coherence times in a room-temperature solid. While the 13 C nuclear 
spin is often a dominant source of magnetic noise, we observe transitions between electron-nuclear 
hyperfine states of NV _ centers in 99% 13 C diamond that are robust to decoherence. At magnetic 
field strengths ranging from 550 - 900 G, these transitions are observable by optically detected 
magnetic resonance (ODMR), and exhibit linewidths narrowed by factors as high as ~130 at room 
temperature over typical electron-type transitions observed from this spin system. We antici¬ 
pate the use of these decoherence-protected transitions, in combination with dynamical decoupling 
methods, for storage of quantum information. 
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INTRODUCTION 


The nitrogen vacancy (NV - ) center is often regarded as a trapped ion in diamond, with 
localized electronic states resulting from its insulating environment and an optical transition 
distinct from the diamond bandgap. Its ground state has a total electronic spin of 1, and 
the spin-spin interaction within the defect results in an energy splitting between the ms = 0 
and ms = ± 1 spin states. Due to optical dynamics and selection rules of NV - centers, the 
ms = 0 spin state fluoresces more intensely than the rris = ± 1 spin states and also becomes 
nearly completely populated upon optical illumination. As a result, the 1115 = 0 spin state 
is referred to as “bright” and the mg = ±1 spin states are referred to as “dark”. This 
spin-dependent fluorescence forms the basis for the method of optically-detected magnetic 
resonance (ODMR) of NV - centers. These properties, in combination with the details of the 
energy level structure of the NV - center, allow for its use in optical sensing of temperature, 
pressure, magnetic and electric fields PHB]; for atomic-scale nuclear magnetic resonance 
(NMR) and magnetic field imaging, as well as for quantum information processing (QIP)[7f- 
H2J. The NV center is competitive as a quantum bit (qubit) candidate due to the ease of 
manipulating its optical and spin degrees of freedom as well as the favorable properties and 
versatile fabrication of the diamond. One significant challenge in the implementation of a 
working quantum computer common to many implementations is producing a qubit that is 
sufficiently decoupled from its environment but may also be interrogated efficiently [13] . In 
this work we address the issue of spin decoherence of NV - centers. 

The sources of decoherence and their interactions with the qubit have been studied for 
many qubit candidates, including photons, phosphorous and silicon nuclear spins in silicon, 
quantum dots, trapped atoms and ions, superconducting circuits, and the NV - center in 
diamond H 2 H5]. Methods of extending coherence times generally include physically re¬ 
moving the sources of decoherence, for example by isotopic enrichment of the host material 
pum, and dynamical decoupling (DD) pulse sequences inspired by the Hahn echo that 
remove interference from static and fluctuating sources of noise [IB]. Another approach to 
achieving useful coherence times is to use transitions that are inherently protected from 
sources of decoherence. One form of these transitions is generally termed Zero First-Order 
Zeeman (ZEFOZ) for satisfying the condition —> 0, in other words when the first-order 
magnetic field ( B ) dependence of a transition with frequency v goes to zero. These are 
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typically found in systems that possess internal interactions (i.e. hyperfine interactions or 
zero field splittings) that do not commute with the Zeeman interaction. For spin qubits with 
such internal structure, magnetic held conditions exist where the Zeeman energy cancels the 
hyperfine or zero-held energy, resulting in an avoided crossing or level anti-crossing (LAC). 
Approaching the regime of the LAC, the magnetic held dependence of spin states becomes 
nonlinear, resulting in ZEFOZ conditions where the transition frequency is desensitized to 
perturbations in the magnetic held. Such phenomena may also be understood by realizing 
that the hrst-order Zeeman shift of a transition frequency is its effective gyromagnetic ratio, 
'Jeff- Spin transitions with weak 7 e ff may be characterized by long coherence times because 
they couple weakly to local magnetic held fluctuations that result in dephasing. This is the 
same reasoning used to explain why nuclear spin coherence times generally exceed electron 
spin coherence times by orders of magnitude: their gyromagnetic ratios are weaker than 
that of the electron by a factor of 10 3 -10 4 . 

The use of ZEFOZ transitions to extend coherence was originally introduced by Longdell 
and coworkers with rare-earth metal ion doped materials [19]. Extended coherence times 
were demonstrated for Pr 3+ :Y2Si05 and Pr 3+ :La2(W04)3 for enhancements of ~630 using 
the ZEFOZ method alone |T9j [20]. Building upon these gains has led to some of the longest 
quantum memory times currently known [21, [22]. Furthermore, the ZEFOZ method is 
broadly applicable. In trapped 9 Be + ions, the ZEFOZ transition between hyperfine levels 
| F = 2, mp = 0) <H-| F = 1, mF = 1) has a coherence time of ~10 s, which is an improvement 
of 5 orders of magnitude over transitions in the same system whose magnetic held dependence 
is linear [23] . Use of a ZEFOZ transition in bismuth-doped natural silicon results in over 
two orders of magnitude increase in coherence time to ~90 ms [23]. Such transitions have 
already been applied to extend coherence of NV - centers in 1.1% 13 C diamond: Lesik and 
coworkers reported an 8-fold increase in coherence time for bulk NV~ centers by tuning 
to a ZEFOZ point at zero magnetic held induced by the hyperfine interaction with 14 N of 
the defect [25] . This gain in coherence was shown to diminish to a factor of 2 for shallow 
implanted NV~ centers or those in nanodiamonds where the defect is subject to stronger 
electric held huctuations. Lastly, Xu et ah employed continuous-wave DD of dressed NV~ 
center spin states, which essentially utilized ZEFOZ transitions in the rotating frame, in 
1.1% 13 C diamond to achieve enhancements of ~20 in spin coherence time [26]. 

In the current study we apply the principles of the ZEFOZ technique to the mixed states 
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of NV - centers in 13 C-enriched diamond and observe a maximum reduction of linewidth by a 
factor of 133.7 from pure electron-type spin transitions typical of the system. We are able to 
explain much of the spectral behavior of this spin system observed by continuous-wave (CW)- 
ODMR using a simple model of four spins: the NV - and its three nearest-neighbor 13 C nuclei. 
Due to our method of detection, we cannot observe transitions directly at ZEFOZ points 
found within this spin system, meaning that while coherence is enhanced significantly, the 
transitions are protected from decoherence (“deocherence-protected transitions” or DPTs) 
rather than true ZEFOZ transitions. 


DECOHERENCE OF NV" CENTERS IN A NUCLEAR SPIN BATH 

The spin dynamics of NV - centers are highly sample dependent, determined largely by 
the spin content of the diamond host as well as by the morphology and treatment of the 
material. Tuning these conditions may result in coherence times as short as hundreds of 
nanoseconds and as long as milliseconds HEII27H1SI]. Consequently, extensive effort has been 
invested in the development of diamond samples ideal for applications that require long spin 
coherence time. 

In bulk single crystalline diamond, the main sources of decoherence for NV - centers may 
be 13 C nuclear spins, PI centers (a spin- 1 /^ electronic defect arising from a substitutional 
nitrogen atom in the diamond lattice), or neutral NV° centers [31]. Because our sample has 
low nitrogen content and high 13 C content, we consider only the effect of the nuclear spin 
bath on NV - linewidth (An). We use this linewidth as an indicator of coherence time (T 2 ) 
by the relation T 2 ~ (An) -1 , where T 2 , the inhomogeneous dephasing time, is less than T 2 . 
Mizuochi et al. reported NV - resonances as narrow as 18 kHz by isotopic purification to 
99.9997% 12 C diamond. In comparison, a dense nuclear spin bath such as that in 99% 13 C 
diamond can lead to linewidths ranging from 50-70 MHz (shown in blue, Figure [lj3) .This 
difference is the consequence of broadening by a range of hyperfine interactions from 13 C 
nuclei occupying all sites surrounding the NV - defect in the diamond lattice. 

We explain some of the behavior exhibited by an ensemble of NV - centers in 13 C-enriched 
diamond (Figure [l}A) in this paper. The ensemble is approximated as a four-spin “molecule” 
composed of the NV“ center and the three nearest-neighbor nuclei, called “first-shell nuclei”. 
The operators S Z1 S x , and S y describe the electron spin angular momentum, individual 
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nuclear spin angular momenta are described by the operators I nx , I ny , and I nz , and the 
total nuclear spin operators are given by K x , K y , and K z . With the applied held along 
the NV - axis, the ground state spin Hamiltonian (H gs ) includes the electron spin Zeeman 
term ('YnvB z S z ), the ground state zero-held splitting interaction ( D gs = 2.87 GHz), and 
the ground state hyperhne Hamiltonian (Hhf^s) of the NV - center interacting with three 
hrst-shell spin-| nuclei, 


Hgs — 7nvB z S z + D 


9S 


S ? 


f C2 
3^ 


+ 11 


HF,gs 


where 

H H F,gs = (S ■ A- h) + RW(S • A ■ ■ A ■ 


. The gyromagnetic ratio of the NV - (7 nv) is ~2.8 MHz/G. We construct the Hamiltonian 
in the frame of the defect axis because our experiment only detects transitions in the NV~ 
spin state along that axis. With the NV“- 13 C hyperhne tensor (H) projected onto the NV~ 
axis, taking the nucleus 13 Ci to lie in the xz plane, the hyperhne term for one nuclear spin 
in HnF,gs becomes 


(S ■ A ■ I\) — A XX S X I\ X + AyySylly + A ZZ S Z I\ Z + A xz (S Z Il X + S x Il z ) 

with hyperhne tensor elements [32j: A xx = 166.9 MHz, A yy = 122.9 MHz, A zz = 90 MHz, 
and A xz = -90 MHz. We obtain the hyperhne interactions with the remaining two nuclei by 
incremental rotations of <p — about the NV~ axis 

The magnetic held dependence of the eigenvalues of the four-spin model is shown in 
Figure [2] This system has a total of 24 eigenstates where each electron spin state is split 
into eight electron-nuclear spin states. We will refer to the sets of eigenstates where the z 
projection of the NV - spin angular momentum is approximately 0, -1, or 1 as the 1115 = 
0, -1, or +1 “manifolds”. In the limit of B = 0, the three hyperhne interactions from the 
hrst-shell spin-| nuclei split the ms = ±1 manifolds each into quartet structures. When 
combining three nuclei of spin- j,, there is only one combination that gives a total nuclear spin 
of ±|, but three possible combinations each that give a total nuclear spin of ±|. This gives 
rise to the 1:3:3:1 intensity ratio of resonances observed in the zero held spectrum, shown 
in Figure [lj3. The transitions leading to this quartet structure (given in blue in the energy- 
level diagram in Figure [lp) are inhomogeneously broadened by the various configurations 
of weak hyperhne interactions of nuclei occupying sites at a distance of 3-8 angstroms 
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[33, M] from the defect. In contrast, the zero held spectrum of an ensemble of NV centers 
in 1% 13 C diamond (shown in black, Figure 03) has one primary resonance, the zero held 
resonance (ZFR), split by crystal strain of approximately 3 MHz. The primary resonance is 
accompanied by two broadened lines of weakened intensity at D gs — 56.9 MHz and D gs + 70.7 
MHz, corresponding to the small probability of finding a NV~ center coupled to a first-shell 
13 C nucleus [35]. It should be noted that the linewidth of this particular ensemble of NV~ 
centers in 1% 13 C diamond is large for NV - centers in diamond of this composition at ~15 
MHz. This is attributed to the proximity of the defects to the surface (no greater than 
100 nm) as surface strain and spin defects are known to produce electric held gradients and 
random fluctuating magnetic helds that lead to rapid dephasing. 

The eigenstates of the nig = 0 manifold are also split into four two-fold degenerate levels 
(see Figure 03) in the limit of B = 0 due to the transverse terms of the hyperhne tensor. 
As a result the structure of the mg = 0 manifold cannot be described by the z-projection of 
total nuclear spin angular momentum. Even at magnetic helds far from any LAC regime, 
these eigenstates exhibit nonlinear magnetic held dependence and sufficient mixing for the 
observation of narrowed DPTs by CW-ODMR spectroscopy. 


RESULTS AND DISCUSSION 

In CW-ODMR of NV~ centers, a transition in NV~ spin state is saturated and observed 
as a decrease in huorescence intensity. The contrast, or normalized difference in fluorescence 
intensity, of these transitions may be explained as a depopulation of the mg = 0 spin 
state [36]. Since mixed eigenstates are not characterized by the pure Zeeman basis, but 
linear combinations of these states, the notion of “bright” and “dark” NV~ spin states 
must therefore be redefined as the spin states with the most or least mg = 0 spin state 
character. Due to small mixing of the electronic spin state in the hyperhne states of the 
ms = 0 manifold, transitions that are nominally forbidden become weakly allowed and have 
sufficient contrast to be detected by ODMR. 

At magnetic held strengths where the influence of electron-nuclear hyperhne interactions 
and the zero-held splitting remains significant, all transitions in the four-spin system are 
possible. In order to gain insight into the nature of the decoherence protected transitions 
observed, we reduce the number of possible transitions to consider by characterizing them 
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with an intensity factor k, which evaluates the probability of their observation by CW- 
ODMR. k is determined by the product of the transition matrix element (TME), a term 
describing the difference in population of final and initial states (A (p)), and a term describing 
the optical contrast (A (S z )) of the transition: 

k = (TME)(A( P ))(A(S 2 z )) 

Where 

TME = {<p f \('y NV (Sx + S y + S z ) + ^(K x + K y + K z ))\ 

A (P> = (<Pf\p\<Pf) ~ {<Pi\p\<Pi) 

A(S'z) = (v f \Sz\<p f ) — (ipi\S z \ipi) 

Here, because optical contrast in ODMR is determined by a change in electron spin state, 
we define the optical contrast term as the difference in expectation value of in the final 
(iff) and initial (< pj ) eigenstates. The transition matrix element is calculated as the inner 
product of the final and initial eigenstates of the spin angular momentum operators with the 
electron (S x , S y , S z ) and total nuclear (K x , K y , K z ) spins, weighted by their gyromagnetic 
ratios. We consider all projections of angular momentum because the orientation of our 
microwave excitation is not exactly known. In this way we do not unintentionally exclude 
transitions that could explain the behavior in the ODMR spectrum. Finally, A (p) is the 
difference in expectation value of the density matrix (p) of the eigenstates involved in the 
transition. To describe the optical pumping, we construct a density matrix where the eight 
ms = 0 sublevels are equally populated and the nuclear spin states are thus unpolarized 
(p = (E — S z )nv <8> Ei3 Cl <g) Ei3 C2 ® Ei3 C3 ). 

In Figure [3j the transitions observed by CW-ODMR spectroscopy are summarized with 
transitions we predict using k greater than 10 ~ 6 , which was determined empirically to best 
fit the data. A more detailed view of the DPTs is given in the inset to the right. Roughly two 
types of transitions may be distinguished: high intensity 1115 = 0 to mg = ±1 electron-type 
transitions (shown in gray), and low intensity DPTs (shown in green). As discussed, transi¬ 
tions of the latter group have much narrower linewidths than the electron-type transitions 
between NV _ spin states due to weakened effective gyromagnetic ratios of the eigenstates 
near LACs, indicating longer dephasing (T 2 *) and potentially coherence (T 2 ) times. 

Two main sets of LACs where DPTs may be found are observed upon inspection of 
the eigenvalues of the NV~-(3) 13 C system in Figure I Set 1 of LACs (Figure |2]C) results 
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from the Zeeman interaction canceling the hyperfine interactions with the first-shell 13 C 
nuclei, and occurs at magnetic field strengths 0-~80 G. Set 2 (Figure IP) results from 
the Zeeman interaction canceling the energy of the spin-spin coupling that leads to the 
ZFS. While a third set of LACs is known to occur from the spin states of the electronic 
excited state near 500 G, the hyperfine tensors of the NV - center and neighboring nuclear 
spins in the electronic excited state are largely unknown. In this study, we develop a simple 
basis for understanding the DPTs that become observable by CW-ODMR at magnetic fields 
approaching but not directly coinciding with the conditions of Set 2 of LACs (~565 - 950 
G), due to the complexity of the behavior of eigenstates in that magnetic field regime (see 
the inset of Figure [3]). The majority of these transitions have energies corresponding to 
transitions between eigenstates in the ms = 0 manifolds. 

It is important to distinguish the transitions we are able to observe in this system us¬ 
ing CW-ODMR from ZEFOZ transitions described thus far in the literature. In previous 
studies, ZEFOZ transitions are observed in spherically symmetric systems only perturbed 
by isotropic hyperfine interactions. For these systems, there exists a “ZEFOZ point” at 
the LAC where the magnetic field dependence of at least one eigenstate involved in the 
transition approaches zero [2Tj. In contrast, our system is perturbed both by hyperfine 
interactions and the spin-spin coupling of the NV“. The transitions we study, which we 
distinguish by terming them decoherence-protected, occur between hyperfine sublevels of 
the ms = 0 manifold, whose magnetic field dependence originates only from mixing with 
ni5 = -1 eigenstates due to the transverse terms of the hyperfine interaction. This means 
that although a ZEFOZ point exists for these transitions at zero magnetic field where the 
coherence properties are optimal, there is a broad range of magnetic fields where transitions 
with enhanced coherence properties may be observed. Furthermore, our detection method 
(CW-ODMR), though highly advantageous for its simplicity, is limited in the sense that it 
requires a transition to have a nonzero change in the z projection of electron spin angular 
momentum (A ms =) for observation. Transitions with the best coherence properties, where 
A ms is closest to 0, are thus not directly observable via optical contrast. 

Nevertheless, we find DPTs between 1119 = 0 sublevels with sufficient optical contrast 
for CW-ODMR at magnetic fields as far from Set 2 of LACs as ~550 G, where A ms 
is on the order of 1CF 3 . The enhancement of coherence is gauged by the linewidths of 
these transitions, which are significantly narrowed in comparison to typical electron-type 


transitions of this spin system. The reduction in linewidth becomes less dramatic as the 
magnetic held approaches Set 2 of LACs, where the 1115 = 0 eigenstates become increasingly 
mixed with those in the nig = -1 manifold and transitions involve increasing A ms- This 
trend is reflected in the nonlinearity of transitions near Set 2 of LACs illustrated in the inset 
of Figure [3j as well as in Figure [4] which shows spectra of DPTs selected over a range of 
magnetic fields (608-871G). It is clear that transitions with increasing optical contrast are 
also more sensitive to magnetic noise. 

The observed transitions are assigned according to energy and the intensity factor k 
of the predicted transition (see Figure [5]). Transitions occurring between the same two 
eigenstates were grouped according to these assignments, and their first- ( 7 e ff = |§) and 
second-order (curvature, C = magnetic held dependence were determined empirically 
using quadratic hts. These parameters for the transitions shown in Figure [4] are given 
in Table [T] and may be used to gauge the transition’s sensitivity to axial magnetic held 
huctuations. In this class of transitions, linewidths as narrow as 527 kHz as well as effective 
gyromagnetic ratios and curvatures as low as 13.62 kHz/G and 0.51 kHz/G 2 are estimated, 
respectively. Comparing with the maximum linewidth of 70.46 MHz for an electron-type 
transition we observe between the ms = 0 and ms = +1 eigenstates away from any set 
of LACs, these transitions are narrowed up to a factor of 133.7. The dephasing times 
estimated from these linewidths may be extended by use of dynamical decoupling methods 

m- 

Upon inspection of Figures [3] and [5j it is clear that the system is too complex to un¬ 
ambiguously assign each transition. Furthermore, direct comparisons between the observed 
linewidths and the effective gyromagnetic ratios characterizing a transition are not always 
consistent. This may indicate homogeneous broadening, which is not taken into account in 
the four-spin model. Inherent in this simplification of the spin system under consideration 
is that all sources of noise causing spectral broadening are axial magnetic fields that are 
static on the time scale of the NV~-(3) 13 C spin dynamics. In reality, the NV“ experiences 
coupling to nuclei occupying all surrounding lattice sites in a volume with a radius as great 
as 5 angstroms from the defect. Experimentally, 15 nuclear sites with interaction strengths 
ranging from 400 kHz to 14 MHz have been measured, whereas 9 nuclear sites with a total 
of 39 symmetrically equivalent positions have been calculated to have interaction strengths 
ranging from 1.5-19.4 MHz [33l 34] . Such a range of coupling provides various relaxation 
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pathways induced by coherent interactions, i.e. nuclear spin flips, or energy conserving 
flip-flops, the understanding of which is beyond the scope of this manuscript. 

Here the general behavior of a class of transitions in a novel spin system has been pre¬ 
dicted. Whereas the amplitude of the transition observed by CW-ODMR becomes stronger 
near Set 2 of LACs, the coherence properties of DPTs approach an optimum far from Set 
2 of LACs, where first- and second-order magnetic held dependence approaches a mini¬ 
mum. These trends are only true, however, for the DPTs we observe using CW-ODMR 
spectroscopy due to Set 2 of LACs. As is evident from Figures [2] and [3j this spin system is 
sufficiently complex to yield DPTs at a variety of magnetic held strengths and orientations. 
Narrow features at low magnetic helds (0.5-600 G, some examples found in the Earth’s held 
CW-ODMR spectrum of NV“ centers in 99% 13 C diamond in Figure [l]) are discussed further 
in a study by Jarmola et. al. |37]. Given the nature of DPTs, where eigenstates are least 
sensitive to inhomogeneous broadening in the limit of (S z ) —> 0, alternate detection schemes 
such as raman heterodyne spectroscopy [38:, [30] may be required to achieve the maximum 
possible extension of coherence time in this spin system. 


CONCLUSION 

We have shown that the complex behavior of nitrogen vacancy centers in 13 C-enriched 
diamond may be approximated by considering the simple spin molecule of the NV _ center 
coupled to the three nearest-neighbor 13 C nuclei. The internal interactions of this spin 
system, i.e. the spin-spin coupling of the NV _ center and hyperfoie interactions with three 
neighboring nuclei result in two regimes of magnetic held where spin states undergo avoided 
crossings and DPTs may be found. Such transitions involve electron-nuclear mixed spin 
states where sufficient projection of NV _ spin onto the mixed state as well as optical pumping 
to the mg = 0 manifold allow for their observation by CW-ODMR spectroscopy. In this study 
we find a subset of DPTs whose linewidths are narrowed as a consequence of transitions in 
spin state being desensitized to the fluctuating magnetic held environment of 13 C-enriched 
diamond. Of these transitions, we estimate effective gyromagnetic ratios as low as 13.62 
kHz/G, with linewidths narrowed by factors as high as 133.7. These results demonstrate a 
method by which NV _ coherence times may be significantly extended, which is a necessary 
step in the development of NV~ centers for use in quantum information processing. 


10 


EXPERIMENTAL 


All spectra were acquired using simple continuous-wave application of laser and microwave 
excitation to detect changes in steady-state populations using a homebuilt confocal mi¬ 
croscopy system. Magnetic fields as high as 1175 G are applied in these measurements by 
mounting a neodymium permanent magnet in a goniometer for control of held orientation. 
The relevant microwave frequencies are applied using waveform generators. The waveforms 
are chopped using a microwave switch and amplified before being sent to a broadband mi¬ 
crowave loop resting directly beneath the diamond sample. Optical excitation is performed 
with a 200mW pumped solid state 532 nm laser. The beam is switched by driving an 
acousto-optic modulator and focused with a 0.7 NA microscope objective to an approximate 
waist of 1.4 /mi (confocal length ~23.6 /.mi, excitation volume ~40 /mi 3 ). After necessary 
filtering and loss due to optics, the power applied to this volume is 4 kW/cm 2 . The fluores¬ 
cence that results is separated from the excitation beam using a dichroic mirror, long-pass 
filtered, and fiber-coupled to a single photon counting module for detection. To describe 
the acquisition of a CW-ODMR spectrum, we refer you to the inset from Figure P- The 
time period for microwave excitation and fluorescence acquisition is 1 ms. The fluorescence 
acquired during microwave excitation is normalized by fluorescence acquired after the NV~ 
ensemble is reinitialized by optical pumping. The sequence for each microwave frequency in 
the spectrum is averaged 2000 times. 

The diamond sample used in these measurements was an electronic grade wafer of di¬ 
mensions 1 x 2 x 0.5 mm grown by Element Six. The sample was grown by chemical vapor 
deposition to a thickness of 500 /mi, where 30 /mi of the thickness on one face was grown 
from 99% 13 C-enriched methane. The estimated concentration of nitrogen impurity in the 
enriched layer is < 1 ppm with less than 10 ppb NV - centers distributed unevenly across 
the sample. These naturally-formed NV _ centers in the 13 C-enriched layer were studied in 
the present experiments. For determination of magnetic held alignment along the NV _ axis, 
an ensemble of nitrogen vacancy centers was created in the opposite face of the diamond 
where the abundance of 13 C nuclei is 1.1%. The creation of this ensemble of NV~ centers 
was accomplished by an implant by Innovion Corporation of 35 keV and 60 keV 14 N ions to 
a huence of 10 13 ions/cm 2 , followed by an anneal under nitrogen atmosphere for 2 hours at 
800° C. The sample is situated over the broadband microwave loop and the focal point of 
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laser beam leaving the microscope objective. A stepper motor with 1 /irn stepping resolution 
in X, Y, and Z dimensions is used to position the implanted ensemble to align the field at 
the focal point of the beam, and then stepped to position the 13 C-enriched layer at the focal 
point for studying DPTs. The quality of field alignment was confirmed by a difference in 
displacement of the UI 5 = 0 to -1 and +1 transitions from the ground state zero field splitting 
of no greater than 5 MHz, allowing us to estimate an error of ± 2 G. The data collected 
at 785 G is an exception, where the difference in displacement was ~10 MHz. In addition, 
we estimate the error associated with using this implanted ensemble of NV - centers as a 
magnetometer by the following equation [ 3 D]: 

5B ~ —- 

7 NV x / N s t m T£ 

Where is the NV - gyromagnetic ratio, N$ is the number of NV - spins, t m is the 
measurement time, and T 2 is the dephasing time. For the concentration and volume of 
excitation given above, N§ — 6.5 x 10 4 . When approximating t m ~ T 2 * ~ (Ay) -1 , where 
(Ay ) -1 is the linewidth of the ensemble given above, we obtain a SB of 0.05G, which is 
negligible in comparison to the field alignment error. 
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TABLE I: Parameters of assigned decoherence protected transitions at 608 G, 739 G, and 871 
G. Center frequencies (n) and linewidths (An) are given in MHz. 7 e ff (^g) and C (f 2 e) are 
given in kHz/G and kHz/G 2 , respectively. For comparison, jnv ~2,800 kHz/G. Observed 
reductions in linewidth (e) are calculated using the average An of electron-type transitions 
observed at magnetic fields away from any set of LACs (65.19 MHz). 
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FIG. 1: The NV-(3) 13 C electron-nuclear spin system. (A) Schematic of the NV - center in 1.1% 
and 99% 13 C diamond. 14 N, 12 C, 13 C, and the lattice vacancy are depicted as blue, gray, dark 
gray, and transparent cyan spheres, respectively. In fully-enriched 13 C diamond, all lattice sites 
around the NV - center are occupied by 13 C nuclei. (B) Earths field ( 0.5G) CW-ODMR spectra 
of ensembles of NV - centers in 1.1% (black) and 99% (blue) 13 C diamond. The inset describes 
the experiment used to acquire all spectra presented in this study (See Experimental). (C) 
Energy level diagram of the NV“ center in the absence of 13 C nuclei in the first shell (black) and 
with three 13 C nuclei in the first shell (blue). The spin quantum number of the NV - is denoted 
by mg, whereas the spin quantum number of the total nuclear spin of the three first-shell 13 C 

nuclei is denoted by m^. 


17 




































FIG. 2: (A) Calculated eigenvalues of the NV _ -(3) 13 C spin system as a function of magnetic 
field strength along the NV - axis. (B) A more detailed view of the 1115 = 0 manifold is shown. 
The transverse terms of the hyperfine tensor result in four doubly-degenerate sublevels at zero 
magnetic field. This degeneracy is broken upon application of a magnetic field, and the eight 
sublevels exhibit a very weak, nonlinear magnetic field dependence from approximately 0-800 G. 

Two sets of LACs exist for the electronic ground state spin Hamiltonian of the NV-(3) 13 C 
system. (C) Set 1 (0.5 - 80 G) occurs between the 1115 = -1 and mg = +1 manifolds and (D) Set 
2 (800 - 1200G) occurs between the mg = -1 and mg = 0 manifolds. 
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FIG. 3: Predicted and observed transitions of the NV - -(3) 13 C spin system. Predicted 
transitions have intensity factors (k) greater than 1CF 6 . This is 67.5 % of all possible transitions. 
The frequency range of the DPTs is magnified at the right for a more detailed view. Many 
weakly-allowed transitions with first-order magnetic field dependence terms (i.e. or 'Yeff 
orders of magnitude lower than qjvy are predicted in this magnetic held regime. 
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FIG. 4: CW-ODMR spectra of DPTs for three different strengths of magnetic held are given to 
demonstrate how the nature of the transitions changes as the magnetic held is increased towards 
the regime of Set 2 of LACs. Gaussian hts of transitions occurring between the same sets of 
eigenstates (W*) are color-coded. Unassigned transitions are not ht. Though the amplitude of 
transitions detectable by CW-ODMR far from Set 2 of LACs is low, the best coherence properties 

are exhibited in this magnetic held regime. 
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FIG. 5: Assigned DPTs (Wj) of the NV~-(3) 13 C spin system compared to all predicted 
transitions at magnetic field strengths near Set 2 of LACs. The observed transitions 
corresponding to the same sets of eigenstates, as well as quadratic fits of their magnetic field 
dependence, are color-coded consistently with those in Figure [4j Quadratic fits are used for 
empirical estimation of first- and second-order magnetic field dependence. The predicted 
transitions corresponding to n > 10~ 6 are underlaid in grayscale. 
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